Non-Gaussianities in New Ekpyrotic Cosmology 
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New Ekpyrotic Cosmology is an alternative scenario of the early universe which relies on a phase 
of slow contraction before the big bang. We calculate the 3-point and 4-point correlation functions 
of primordial density perturbations and find a generically large non-Gaussian signal, just below the 
current sensitivity level of CMB experiments. This is in contrast with slow-roll inflation, which 
predicts negligible non-Gaussianity. The model is also distinguishable from alternative inflationary 
scenarios that can yield large non-Gaussianity, such as DBI inflation and the simplest curvaton- 
like models, through the shape dependence of the correlation functions. Non-Gaussianity therefore 
provides a distinguishing and testable prediction of New Ekpyrotic Cosmology. 



The ekpyrotic scenario is a candidate theory of early- 
universe cosmology. Instead of invoking a short burst of 
accelerated expansion from a hot initial state, as in infla- 
tion, the ekpyrotic scenario relies on a cold beginning fol- 
lowed by a phase of very slow contraction. Despite such 
diametrically opposite dynamics, both models predict a 
flat, homogeneous and isotropic universe, endowed with a 
nearly scale invariant spectrum of density perturbations, 
and are, therefore, equally successful at accounting for 
all current cosmological observations. 

An important drawback of the original ekpyrotic the- 
ory [l| is how to avoid the big crunch singularity without 
introducing ghosts or other pathologies. Moreover, the 
fate of perturbations through the bounce is ambiguous 
— whereas the scalar field fluctuations are scale invari- 
ant during the contracting phase, the curvature pertur- 
bation on uniform-density hypersurfaces, £, is not. And 
since the latter remains constant on super-horizon scales, 
one generally expects its (unacceptably blue) spectrum 
to be preserved irrespective of the bounce physics. De- 
spite considerable work indicating that stringy effects at 
the bounce may positively alter this conclusion (2j, the 
issue of matching conditions remains controversial. 

Both of these issues have been resolved in the re- 
cently proposed New Ekpyrotic scenario !■ In 1, we 
derived a fully non-singular bounce within a controlled 
and ghost-free four-dimensional effective theory using the 
ghost condensation mechanism Moreover, £ acquires 
a scale invariant spectrum well before the bounce, thanks 
to an entropy perturbation generated by a second scalar 
field 0, S, Hf. Thus New Ekpyrotic Cosmolo gy appears 
to be a consistent alternative to the inflationary scenario. 

A distinguishing prediction lies in the tensor spec- 
trum [lj: inflation predicts scale invariant primordial 
gravity waves, whereas ekpyrosis does not. Detecting 
tensor modes from CMB B-mode polarization could rule 
out the ekpyrotic scenario, whereas an absence of detec- 
tion would not discriminate between the two models. 

In this Letter we focus on another key observable: the 
non-Gaussianity of primordial density perturbations. We 
show that New Ekpyrotic Cosmology generically predicts 
a large level of non-Gaussianity, potentially just below 
current sensitivity levels and detectable by near-future 



CMB experiments. 

We calculate the 3-point and 4-point functions. For 
typical parameter values, the amplitude of the 3-point 
function is generically large, with /nl around the current 
WMAP bound Q: -36 < / NL < 100. That is, assuming 
all parameters are 0(1), /nl approaches the limits of this 
bound, depending on the sign of a parameter. These val- 
ues are well above the expected sensitivity of the Planck 
experiment: |/nl| 20. The amplitude of the 4-point 
function is also generically large: tnl ~ 10 4 , which is 
again near the estimated WMAP bound and within the 
reach of Planck: r NL < 600 @|. 

This is in stark contrast with the highly Gaussian spec- 
trum predicted by slow-roll inflation. Comparably large 
non-Gaussianity does arise in non-slow roll models, such 
as DBI inflation Q , and whenever the precursor of den- 
sity fluctuations is a light spectator field, such as in the 
curvaton [l(| H3 or modulon scenarios [IH, [l^|. How- 
ever, as we will see, New Ekpyrosis predicts a different 
shape dependence in momentum space for the 3- and/or 
4-point spectrum than the simplest such models. 

Non-Gaussianity therefore offers a distinguishing pre- 
diction of New Ekpyrotic Cosmology, potentially testable 
in CMB experiments within the next few years. 



I. NEW EKPYROTIC COSMOLOGY 

As with inflation, ekpyrosis relies on a scalar field <f> 
rolling down a potential V (</)). Instead of being flat and 
positive, however, here V((f>) must be steep, negative and 
nearly exponential in form. For concreteness, we take 



VW>) = -V e'^ A 



(1) 



where A = y/eMp\ and K<1. The Friedmann and scalar 
field equations then yield a background scaling solution, 



a(t)~(-t) 2 *; 4>(t) = A log 
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with Hubble parameter H = 2e/t. Since e -C 1, this 
describes a slowly-contracting universe with rapidly in- 
creasing H, again in contrast with the rapid expansion 
and nearly constant H in inflation. 
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In single-field ekpyrosis, fluctuations in 4> acquire a 
scale invariant spectrum. As we review shortly, this 
traces back to the fact that the above solution satis- 
fies V ; 00 = —2/t 2 . However this contribution exactly 
projects out of £, leaving the latter with an unacceptably 
blue spectrum. Since £ is conserved on super-horizon 
scales barring entropy perturbations, it is generally ex- 
pected to match continuously through the bounce, al- 
though stringy effects could alter this picture Q. 

New Ekpyrotic Cosmology introduces a second field, x, 
as the progenitor of the scale-invariant perturbation spec- 
trum 1 3, 5]. This field has no dynamics during the ekpy- 
rotic phase and remains approximately fixed at \ — 0. 
However, as we describe below, its fluctuations generate a 
scale-invariant spectrum of entropy perturbations, which 
gets imprinted onto C at the end of the ekpyrotic phase. 

An essential condition in obtaining a scale-invariant 
spectrum is that at x = the curvature of the potential 
be nearly the same along the x an d 4> directions: V^ 
V^. An example of such a potential is 

n^x)=v(0)(i+^+3 fS 3 + ¥XI 

(3) 

The higher-order x terms are naturally expected to be 
suppressed by the same scale A as the quadratic term, 
hence the form ([3|). For simplicity we take a.%, an, . . . to 
be constants. While potential © yields a slightly blue 
spectral tilt, a more general potential is presented in Q 
which allows for the observed red tilt without altering the 
conclusions for non-Gaussianity arrived at in this paper. 
Note that the required field trajectory lies along an un- 
stable point. However, a pre-ekpyrotic, stabilizing phase 
can easily create initial conditions so that this trajectory 
is arbitrarily close to the tachyonic ridge 

Power spectrum for x' Since our space-time back- 
ground is nearly static, we ignore gravity in studying x 
perturbations. To linear order, the Fourier modes Sx^ 
around x = satisfy a free field equation with time- 
dependent mass V xx = V<p<p = —2/t 2 : 



Sxk 



(0) 



Assuming the usual adiabatic vacuum, we find 

— ikt 
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On super-Hubble scales, k(—t) <C 1, the power spectrum, 
defined by (Sxk^Xk?) = (^) 3 S a (k + k')P x (k), is 



3 1 



(6) 



which is scale invariant. Including gravity and departing 
from the pure exponential form (fT]) results in small devi- 
ations from scale invariance. This can yield a small red 
tilt, consistent with current CMB observations [3j. 



Evolution of £: We focus for simplicity on the regime 
where all relevant modes are well-outside the horizon, 
k <C aH. In the small-gradient approximation, the met- 
ric can be written as ds 2 = -Af 2 dt 2 +e 2 ^ S ' t " l a 2 (t)dx 2 fl4| . 
where Af is the lapse function, and £ is the curvature 
perturbation. The evolution of £ on uniform-density hy- 
persurfaces is governed by 



( = 2H- 



SV 



25V 



(J) 



where SV = V(4>, x) ~ V(<f>,x)- A key simplification is 
that 8(j> has a steep blue spectrum at long wavelengths 
and, hence, can be neglected. Thus, for the potential ([3]), 
we have SV « V(4>)S X 2 /2A 2 + ... 

To proceed further, one needs an expression for Sx to 
higher-order than the "free" part Sx^ ■ To do this, we 
solve Sx + Vxx^X — 0, valid at long wavelengths, pertur- 
batively: Sx — Sx^ + Sx^ + ■ ■ ■ To lowest order, this 
equation reduces to ^ in the limit k — > 0. The next or- 
der, S x {1 \ satisfies S x {1) + V XX 5 X (1) + V XXX (S X {0) ) 2 /2 = 
0. Using ©, © and Sx {0) ~ l/t, we find 



sx=5x^+^(s x ^y+... 

Substituting into ([7)1. one can integrate to obtain 



Cok - 2 
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The ekpyrotic phase must eventually end if the uni- 
verse is to undergo a smooth bounce and reheat into a 
hot big bang phase. This is achieved by adding a feature 
to the potential (|3|) which eventually pushes x away from 
the tachyonic ridge Q. Denote the time at which ekpy- 
rosis stops as t en d- For simplicity, we model this with 
V x suddenly becoming non-zero and nearly constant at 
X = 0. Denote this constant by V x \. The exit phase 
is assumed to last for a time interval At which is short 
compared to a Hubble time: \H en< x\At <C I . This will be 
the case provided the potential satisfies 



v x \ 2 ~ • 



(10) 



The exit phase generates an additional contribution to 
£. To compute this in the rapid-exit approximation, we 
can treat the right-hand side of {2} as approximately con- 
stant. In evaluating this constant, note that, to leading 
order in Sx, we have SV « U x |c>x = ±H 2 nd MpiSx/ \[£ 



(Higher-order terms in Sx yield small corrections to ([9]) 
and are therefore negligible.) Thus £ changes from £ e k 
by an amount £ c during the exit, given by 



SxiUnd) 



pi 



(11) 



where (3 = \H en ^\At-\/e/e x . Noting that to lowest or- 
der Sx ~ Sx^ ' an d substituting I© evaluated at i on d, it 



3 



follows from (fTTj) that the £ power spectrum is 



k 3 P c (k) 



Ae(3 2 



k 3 p x {k) = p 
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Up to the prefactor /3 2 , this is identical to the inflation- 
ary result, with e playing the role of the usual slow-roll 
parameter. In the exit mechanism of (3 denotes the 
overall change in angle in the field trajectory: (3 = A9. 

Let us pause to discuss the parameter values that sat- 
isfy the CMB constraint k 3 P ( (k) w lCT 10 . Although H 
passes through zero at the bounce, as argued in [3J] its 
magnitude is essentially the same at the beginning of the 
hot big bang phase as it was at i on( j, the end of the ekpy- 
rotic phase. In other words, -ff e nd sets the reheat tem- 
perature in the expanding phase. For GUT-scale reheat 
temperature, we have H cn a/Mp\ w 10~ 6 . Meanwhile, 
(3 is a free parameter whose value depends on the exit 
dynamics. For the explicit exit mechanism of 3], how- 
ever, the natural value is (3 ~ 0(1). In this case, setting 
k 3 P^(k) = 1CU 10 implies e w 10~ 2 . We will henceforth 
take (3=1 and e = 10~ 2 as fiducial parameter values. 

Combining (JTTJ) with © and © yields 



(( x )=( c ( x ) + —C 2 c (x)T 



5a 3 



U4e 2 (3 3 



C c 3 (x) + 



(13) 



The exit from the ekpyrotic phase is followed by a ghost 
condensate phase which leads to a non-singular bounce 
and reheating. Meanwhile, x S e ^ s stabilized and further 
evolution is governed by the single scalar <p. It follows 
that C is conserved through the bounce and emerges un- 
scathed in the hot big bang phase. 



II. NON-GAUSSIANITY 

3-point function: The 3-point £ correlation function in 
New Ekpyrotic Cosmology is given by Q 

(Ck i a 2 Ck 3 ) = (27r) 3 S 3 (k 1 +k 2 + k 3 )B(k 1 ,k 2 ,k 3 ), (14) 
where the shape function B{k\, k 2 , k 3 ) is 

B(k u k 2 , k 3 ) = ^/ NL {P c (fc 1 )P c (fc 2 ) + perm.} . (15) 

This is of the so-called local form [l5j]. Equations (fl"4"l) 
and (fT5"f are consistent with ((x) of the form ((x) = 
Cg(x) + ^fNLCg( x ), where ( g is Gaussian. The correlation 
function is evaluated at £ en( j ignoring gravity. 

Thus the 3-point function is fully specified by /nl ^16]. 
This parameter receives two contributions. To begin 
with, the non-Gaussianity of <5%, due to its cubic inter- 
action in ([3]), is inherited by ( through (fTTj) . Following 
Maldacena [18[ , the 5x 3-point function is given by 

(5X1^X2^X3) = -i ds(0\ [6x^X2^X3, H in t(f 



(16) 



where 8xi = $x( x i)-> an d "Hint is the cubic interaction 
Hamiltonian from ©: H int = V(0)a 3 x 3 /3!A 3 . An ex- 
plicit calculation yields the intrinsic contribution 



j?int 
J NL 



5 a 3 
24 (3t 



(17) 



The =p sign corresponds to choosing V x \ to be ±. 

The second contribution comes from the non-linear re- 
lation between 6x and C embodied in (jTTJ) and (TT3]). Even 
if Sx were Gaussian, this non-linearity would make £ non- 
Gaussian. This conversion contribution to /nl is: 



J'COllV 

/NL " 24 (3 2 e • 
Summing (fT7|) and (fT8|) yields a combined f\ 



NL- 



r f'mt I rcorc 

/NL = /NL "+" /NL 



24/3 2 e 



(1 T a 3 0) 



(18) 



(19) 



Since this is inversely proportional to e <C 1, non- 
Gaussianity tends to be large in New Ekpyrotic Cosmol- 
ogy. Related ekpyrotic models @,[l7[ also give /nl ~ 
(A ghost condensate bounce and second scalar field are 
also invoked in Q, albeit without an explicit conver- 
sion mechanism; and while the two-field ekpyrotic phase 
of (l7| is similar to ours, the bounce physics remains un- 
specified.) This is in sharp contrast with slow- roll infla- 
tion, where /nl is proportional to the slow-roll param- 
eters and therefore unobservably small. For concrete- 
ness, consider our fiducial model with GUT-scale reheat- 
ing, (3 = 1 and e = 1CP 2 . Taking, for example, the — 
sign in (fT?) and choosing 2.728 > a 3 > —3.8 yields /nl 
within the present WMAP 2a range: -36 < /nl < 100. 
Thus a 3 ~ 0(1) yields a non-Gaussian signal near the 
WMAP bound. Lower reheating temperatures corre- 
spond to smaller e and, therefore, larger non-Gaussian 
signal. Of course, |/nl| can always be made smaller by 
taking (3, e to be larger and/or by suitably choosing a 3 . 
4-point function: The connected 4-point function, 

(CfciCtaCfcaCfc) = (2n) 3 S 3 (k 1 +k 2 + k 3 + k 4 ) 

■ [Ti^ , k 2 , k 3l ki) + T'(fci , k 2 , k 3 , k A )}, (20) 

involves two different shape functions, evaluated at £ e nd : 

T = irNL{n(fci)P C (fc 2 )P c (fci 4 ) + 23perm.} ; 
T' = K NL {P c (k 1 )P c {k 2 )P c {k 3 ) + 3 perm.}, (21) 



where fc, 



Thus T and T" are specified re- 



spectively by the tnl and knl parameters. (Note that 
knl is proportional to the f 2 parameter of Q.) Equa- 
tions and (|2ip are consistent with ((x) of the form 



where („ is Gaus- 



sian. Note that we can obtain tnl immediately by simply 
comparing its definition with that of /nl : 



36 f2 
25 /nl 



1 



16/3V 



(lTa 3 (3f 



(22) 
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This was also checked by explicitly computing the three 
and four-point functions. 

Let us now consider knl- It receives two contribu- 
tions: i) an intrinsic piece due to cubic and quartic §x 
interactions; ii) a conversion piece from the non- linear 
relation between 8\ and C- The first contribution arises 
from cubic and quartic terms in \ in the potential ([3]). 
An explicit calculation gives 



K- mt 



2a 4 + 3a 2 . 
40/? 2 e 2 



(23) 



The second contribution is encoded in the £ 2 and (j^ 
terms in (fT3"|). Comparing with (|2H)) . we obtain 



5a 3 



"NL 



1 24/3 3 e 2 ' 

Combining the above results, we find 



„int 



"NL 



a 3 (9a3/3T25) + 6Q4/? 
120/3 3 e 2 



(24) 



(25) 



Both tnl and are inversely proportional to e 2 and 
therefore also tend to be relatively large. Note that tnl 
is always positive, whereas knl can be positive, zero or 
negative depending on the choices of a 3 and a^. For 
instance, our fiducial parameter values for GUT-scale 
reheating with a^^a^ ~ 0(1) yield tnl ~ 10 4 , which 
is around the estimated bound for the WMAP experi- 
ment |8j. Lower non-Gaussianity can again be achieved 
by taking larger (3, e and/or by a suitable choice of a 3 
and a.4. 

Discussion: The simplest inflationary models, consist- 
ing of one or more slowly-rolling scalar fields, all predict 
negligible 3-point and higher-order correlation functions. 
Non-Gaussianity therefore offers a robust test to distin- 
guish New Ekpyrotic Cosmology from slow-roll inflation. 



Significant inflationary non-Gaussianity can be ob- 
tained in non-slow- roll models, such as DBI inflation, 
albeit with a distinguishable shape dependence. Our 3- 
point function is "local" , characterized by a momentum 
dependence that peaks for squeezed triangles, whereas 
the DBI amplitude peaks for equilateral triangles [l5| • 

Large non-Gaussianity may also be achieved in the cur- 
vaton scenario. The curvaton 3-point function is also 
of the local form and hence cannot be used to distin- 
guish curvatons from New Ekpyrotic Cosmology. There 
is, however, an essential difference at the 4-point level. 
In the simplest curvaton model, the progenitor of density 
perturbations is a free field. Thus, knl ~ /nl [HI- In 
contrast, in New Ekpyrosis, tnl and knl are generically 
of comparable magnitude (~ e -2 ) and are expected to 
exhibit a distinguishable shape dependence. More intri- 
cate curvaton models with self-interactions can also yiel d 
large «nl- Similarly for general modulon scenarios [13| . 

Near-future non-Gaussianity observations will, there- 
fore, test the new ekpyrotic paradigm and can potentially 
distinguish it from its inflationary alternatives. 

In this paper we have used a simplifying approxi- 
mation so as to obtain an analytic expression for non- 
Gaussianity. Our results give the exact parametric de- 
pendence while the details of the potential, roll-off time 
and so on are encoded in model dependent parameters, 
such as p. We have checked that different quasi-analytic 
approximations continue to give the same parametric 
dependence as presented here, although with model- 
dependent coefficients that can at most differ from those 
in this paper by factors of order unity. 
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